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Abstract 

Biological responses often obey Weber's law, according to which the magnitude of the response 
depends only on the fold change in the external input. In this study, we demonstrate that a system 
involving a simple autocatalytic reaction shows such response when a chemical is slowly synthesized 
by the reaction from a faster influx process. We also show that an autocatalytic reaction process 
occurring in series or in parallel can obey Weber's law with an oscillatory adaptive response. 
Considering the simplicity and ubiquity of the autocatalytic process, our proposed mechanism is 
thought to be commonly observed in biological reactions. 

PACS numbers: 82.39.-k, 87.10.-e, 05.45.-a 
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Adaptation is ubiquitous in biological systems. One kind of adaptation is perfect adap- 
tation, wherein some state variables of a biological system change in response to a change 
in the external conditions but slowly come back to their original values Such adaptation 
is observed in various systems such as signal transduction systems, sensory systems, and 
neural systems. 

Several theoretical studies have investigated such adaptation by using a set of ordinary 
differential equations, in which certain variables return to their original value independently 
of the parameter values representing the external condition. Several models for studying 
perfect adaptation have been proposed and analyzed 2H8J. 

Responses in a biological system have another ubiquitous characteristic in addition to 
adaptation: the degree of response is often proportional to the ratio of an external stimulus 
before and after it is applied, rather than the absolute difference. In other words, the 
response detects the fold change, which is commonly known as Weber's law. This kind of 
response was discovered in the field of psychology, wherein it was observed that our sensory 
response is proportional to the logarithm of the magnitude of an external stimulus. 

Recently, Weber's law was also observed to be obeyed by cellular responses in microor- 



ganisms 



9l-|ll|. Shoval et al. [l2| further defined fold-change detection (FCD) in a sense 
stronger than Weber's law: not only the magnitude of the response peak but also the entire 
relaxation profile over time depends only on the change ratio of the stimuli. In other words, 
the relaxation profile over time is identical as long as the ratio is constant, irrespective of 
the absolute magnitude of stimuli. Shoval et al. also developed a theoretical model for 
such FCD by using a feed-forward gene regulation network. However, it is currently unclear 
whether this FCD is ubiquitous, whereas Weber's law itself seems to be quite universal in 
biological responses. 

An adaptive response usually has two components: rapid response to a change and slow 
relaxation to the original state. Oosawa and Nakaoka confirmed the existence of these 
distinct timescales in the adaptive response of Paramecium to chemotaxis and also the 
relevance of their existence to chemotaxis [3] (see also {?]]). It is therefore important to 
elucidate the relevance of such a timescale difference to Weber's law for adaptation. 

Most biochemical processes involve catalytic reactions, whereas the growth and repro- 



duction of a cell involve an autocatalytic process 



16| . In the present Letter, we show 



that a simple system with catalytic reactions exhibits an adaptive response according to 
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Weber's law if the autocatalytic process is slow [3, H]- Because the autocatalytic process 
often involves several steps and requires a long time for completion, the ubiquity of adaptive 
responses obeying Weber's law of a cell can be easily understood from this simple system. 
We also demonstrate Weber's law in a system of chained or parallel-connected autocatalytic 
reactions. In addition, we explore oscillatory adaptation with Weber's law and establish a 
condition for it. Finally, we discuss the relevance of the results to biological responses. 

First, we study a simple reaction model of two chemicals Xq and X\ with the following 
autocatalytic reaction, as introduced in X + X\ — > 2X\\ the model also includes (a) the 
synthesis of X from an external resource chemical S as S — > X and (b) the degradation 
of Xq and X\. By representing the concentrations of the two chemicals as Xq and x\ and 
suitably scaling the time and concentration variables, we get the rate equation as 

dxo dx\ xqXi — x\ ... 

= S-x x 1 -x , — = . (1) 

at dt r 

The steady state is given by = S— 1. According to a linear stability analysis, this 

state is stable when S > 1, i.e., as long as x\ > 0. It should be noted that Xq is independent 
of S. The chemical concentration responds to the concentration of the external signal, S; 
when S increases (decreases), xo increases (decreases) from its steady-state value but later 
returns to the original value. Thus, X always shows an adaptive response, whereas the 
steady-state concentration of X 1 changes with changing S. 

To examine Weber's law, we study the response of xq when S changes as So —> pSo(p > 0), 
and calculate the dependence of the peak value of x during this change in S. We assume 
r > 1, which is required to ensure a fast response and slow adaptation. Under the adiabatic 
limit, x\ changes more slowly than xq does. Then, during the fast response of xq to the 
change in S, X\ can be assumed to remain at the steady-state value under the condition 
of S — So- Then, the peak value of x is obtained from (dx /dt) _ pea k = by fixing the 

' Xq — Xq 

value of X\ to x\ = So — 1. A straightforward calculation gives us x v § ak = p. Hence, xq 
changes from the original value Xq — 1 to x P Q ak = p and then returns to the original value. 
This amplitude of the response depends only on p, i.e., the ratio of the shift in S, and is 
independent of the So value. Thus, our model demonstrates Weber's law under the adiabatic 
condition. 

From a standard linear stability analysis, we get two timescales for this adaptive response: 
one for response and the other for adaptation. When S is large enough, we can assume that 
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FIG. 1: Plot of ratio of maximum response amplitude Ai/Aq (ordinate axis) versus r(abscissa 
axis). Aj = x^ akt — Xq with the change Si — > pSi. Ai/Ao = 1 implies the independence of the 
response amplitude from the Si value, or Weber's law. Data with three different (So, S\,p) sets are 
shown: (S ,Si,p) = (5,50,2) for +, (S ,S x ,p) = (5,500,2) for x, and (S ,Sx,p) = (5,50,10) for 
*. 



the timescale for adaptation is given just by ~ r and is independent of S. However, the 
timescale for response (the peak time) is given by ~ 1/5*, which is still dependent on S. 
Thus, the temporal profile of the response depends not only on p but also on So. In this 
sense, our considered model does not satisfy the condition for FCD as defined by Shoval et 

ai. IrJ. 

We note that the adiabatic condition need not be rigid in order for Weber's law to 
be demonstrated reasonably well. Fig JTJ shows a plot of the peak- value ratio x^ eakl /xQ eak2 
obtained by multiplying S p— fold from Sq or Sq, respectively, as a function of r. The 
ratio is close to unity; it is independent of the initial So as long as r > 10. In a previous 



study, we determined that for efficient chemotaxis 



a, 



the ratio of the response time to sense 



n 

microorganisms such as Paramecium [13] and E. coli 



Tn) 
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19] . In the present model, the 



timescales for response and adaptation are given by r s ~ 1/S and r a ~ r, respectively, and 
the above condition corresponds to r ~ 10. For this timescale ratio, the peak-value ratio is 
~ 1.2 as shown in FigJU and thus, Weber's law is approximately obeyed. 

Adaptation obeying Weber's law is not restricted to this two-component reaction system; 
in fact, it is generally observed in a system with autocatalytic reactions with a slower 
timescale. Here, we study several extensions of autocatalytic reaction networks that show 



adaptation obeying Weber's law. First, we consider an autocatalytic reaction occurring in 
series with N chemicals, as S — >■ X — y X\ — > ... — > X^-i- The rate equations are given by 

dxo/dt=( S — P\XqX\ — xo)/tq, 
dxi/dt = (fcxi-iXi - f3 i+1 XiX i+1 - Xi)/Ti (i^0,N - 1), 
dx N -i/dt = ((3n-iXn-2Xn-i - x N -i)/tn-i- (2) 

Then, the steady-state concentrations are obtained as Xq = S/((3\x\ + 1), x* = (fii+2%* + 2 + 
(i 7^ 0, iV — 2), and x* N _ 2 = 1//3n-i- Accordingly, both the steady-state con- 
centrations x* N _ 2 and x* N _ 2m (m > 1) are independent of the external-signal concentra- 
tion S, and these results show the adaptive response to the change in S. Here, we also 
note that the adaptation alone of the chemical x^-2 depends only on the reaction process 
Xns X N ^ 2 -Xjv-i- As long as this chain reaction is autocatalytic and X N _ 2 increases 
with S, xn-2 undergoes perfect adaptation. We can modify other reaction processes while 
retaining the adaptive response of Xn-2- 

First, we study the case of = 3 in detail. For simplicity, we use = 1. Suppose Xi > 0; 
the steady state is given by Xq = S/2, x\ = 1, and x\ = (S — 2)/2, which is linearly stable 
if S > 2 and X\ shows adaptive response to the change in S. 

The relaxation process to the state, however, depends on the timescale, as is con- 
firmed from the eigenvalues of the Jacobi matrix (FigfS]). To achieve a normal adaptive 
response with monotonic relaxation, Tp -C Tl,t 2 is required (for large S, 8r < r 2 and 



To/2 < ti/S < r 2 /32 are required [20(]): otherwise, the system would show an adaptive 
response with damped oscillation. 

An adaptive response with damped oscillation has often been observed in experiments, 



e.g., in [21 



22]. In contrast, theoretical models for such oscillatory adaptation have rarely 
been explored. Our proposed model provides a simple example of an oscillatory adaptive 
response. 

In an autocatalytic reaction chain as well, Weber's law for responses is obeyed for only 
a certain range of parameters. Here, we study the conditions for Weber's law in the N = 3 
case as an example. To analyze this case, we again study the system response to the change 
in S as So — > pSo(p > 0). By setting r < ri C r 2 , we study the dynamics of X\ and 
its peak value (xf eafc ) during an adaptive response. Since the change in x is much faster 
than that in x±, we can assume that Xq always takes the equilibrium value defined from 
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FIG. 2: Behaviors of an adaptive variable(xi) in the N = 3 case. Responses corresponding 
to S = 100 ->■ 200 at t = and S = 200 -> 400 at i = 40 are shown respectively. (Left) 
tq = 0.01, ti = 1,T2 = 10 and (Right) tq = 2, n = 1,T2 = 1, zoomed in the inset. 



the present x\ value, and hence, x is obtained from dx /dt = as x = pS /(l + X\) for 
given x\. However, X2 still remains at the equilibrium value under the condition S = So, 
i.e., X2 = (So — 2)/2. Using these values, x^ eafc that satisfies (dxi/dt) xi=3 j>eak = is given by 
x± eak — 2p — 1. Accordingly, x^ eafc — = 2(p — 1) depends only on the ratio of change in S 
and not on the So value itself. Thus, the model given by eq.(j2J) with N = 3 obeys Weber's 
law under the limits of t /ti — >■ and t 2 /ti — > oo. 

The above condition for Weber's law is a straightforward extension of the condition 
derived in the original model with two variables. To examine the validity of the above 
approximation, we numerically studied the peak-amplitude ratio of x\ under varying t 2 and 
t and constant T\. This condition for all r's corresponds to the upper-left area in Figj3j 
where Weber's law is confirmed. In this region, x\ shows a monotonous relaxation. 

In Figj3j however, we find another regime that (approximately) satisfies Weber's law 
in the lower-right triangle, where x\ shows a damped oscillation and the above adiabatic 
conditions are not satisfied. Indeed, in this case, Weber's law is well obeyed for large S and 

We now consider the response in the case that S changes as So — > pS (p > 0) with 
t 3> Ti,t 2 and S 3> 1. We solve eq.© with the initial condition x (0) = So/2, Xi(0) = 1, 
and #2(0) = So/2 — 1 to obtain the first peak in X\ at t — t peak , which is expected to be 
t peak <c 1 because S is sufficiently large. From the assumption of large S and To ^> T\, xq is 



6 



10000 
100 



C\J 



1 



0.01 
0.0001 



I I I I 


"'I I I I I ! 


' I I ... 





1.4 
1.3 
1.2 
1.1 
1 



0.0001 0.01 1 100 10000 

T0/T1 



FIG. 3: Phase diagram with tq/t\ (abscissa axis) and t 2 /ti (ordinate axis). Ratio of the maximum 
response amplitude, A1/A0, is shown. Ai = x^ ak — x* with S = 200 — > 400 and Ao with 
S = 100 ->• 200. 



solved approximately as xo(t) ~ ^ + ^ T * ' S ° t for small t. Then, we get 

dlnxx S Q (p-l) S , , , . , , 

—r. 1-— — {x 2 (t -x 2 }. 3 

at t tx (S - 2)ri 

Temporarily neglecting the latter term for small t, we get Xx ~ exp [S (p — l)t 2 /(2r r 1 )]. 
By substituting this in d(\nx2)/dt = (So — 2)(xx — l)/(2r2), we get 22 = 
£2(0) exp [(p — 1) Sot 3 1 (6T0T1T2)]. Now, with this increase in £2, the second term in eq.(j3]) is 
no longer negligible, as a result of which the sign of eq.(|3]) changes to negative and Xx stops 
increasing. Thus, the time of occurrence of the first peak in xx is estimated from the time 
when the first term is equal to the second one as 

{tpeak) 2 „ i 2ti t 2 /So. (4) 

By substituting this value in the expression of Xx, the peak value of Xx is roughly estimated 
as 

6(p - l)r 2 



peak 

x\ ~ exp 



(5) 

r o J 

Hence, the maximum amplitude depends on the change ratio of S and not on Sq. Next, we 
numerically verified the dependence of the peak amplitude on r's from the original equation 
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FIG. 4: Comparison between results of numerical simulation (points) and approximate calculation 
(lines). ln(xf afc ) and (fP eak ) 2 are shown as a function of t 2 with t\ fixed to 1. ln(X] )eafc ) with 
r = 100(+) and r = 1000(x) and (tP eak ) 2 with r = 100(*) and r = 1000(D) are taken from the 
numerical simulation. Eq.([5]) with tq = 100 (dotted line) and tq = 1000 (broken line) and eq.(j3]) 
(straight line) are also shown. 



(eq.(j5])). As shown in FigHJ the above approximations give good estimates for the peak 
position and explain the validity of Weber's law. 

Note that we need to satisfy t peak <C r as r takes the largest value among the three in 
this area. Then, we get 

TlT 2 S . . 

If <<: 12' (6) 

which explains the boundary condition in the lower-right area. Note that the peak value 
itself increases as t 2 /tq increases. Hence, to achieve an oscillatory response with a large 
amplitude, the value on the left-hand side of eq.(ll]) should preferably be large. Then, an 
oscillatory response with Weber's law with a relatively large peak will be obtained at the 
border of the above condition. 

For iV > 3 as well, xn-2 shows an adaptive response with both monotonic relaxation 
and damped oscillation. For the former case, by carrying out similar analysis with N = 3 
case, Weber's law is shown to hold under the condition of r , • • • , t^v— 3 tn-2 "C Tzv-i- 
According to eq.(j2]), XN-2m {m > 1) also shows an adaptive response, and Weber's law for 
%N-2m is satisfied when r , ■ ■ • , r N _ 2m -i "C r N _ 2 m T N-2m+i, i-e., if all processes of influx 



S 



are fast and only the next runoff process is slow. Weber's law with damped oscillation is 
also satisfied for large S and r . 

We also studied the case of catalytic reaction networks having parallel paths. Two types 
of extensions are possible. In one extension, many terminal chemicals exist instead of a 
single X\ in eq.flTJ, which are reached in parallel from xq (S — > Xo — > x\, x\, x\, ■ ■ ■ x\). It can 
be straightforwardly shown that variable xq shows an adaptive response, since the mean field 
of all x{'s acts as the variable of x\ in the original two- variable case. Moreover, in this type 
of extension, each catalytic reaction need not be autocatalytic; that is, the reaction from xq 
to x{ can be catalyzed by x™ with j ^ m. Adaptation obeying Weber's law is possible in 
the adiabatic limit. The other extension includes reaction paths in two or more rows, each 
of which is a series reaction, as in the case of variable iV (S — > y — > a\ — > ■ ■ ■ a N — y y 1 and 
S — y Do — y b± — y ■ ■ ■ bj\j — y yi). Here, both an input variable (yo) and an output variable 
(yi) can show adaptive responses only when there are an odd number of elements in each 
row (N and M). In this case, adaptation obeying Weber's law is again possible, but there 
usually are more restrictions on the parameter values. 

In this Letter, we have demonstrated that a simple autocatalytic reaction process in a 
system leads to its adaptation that obeys Weber's law. We first confirmed such an adaptation 
in an autocatalytic reaction of two variables. Coupling of one more variable to the reaction 
led to oscillatory adaptation obeying Weber's law. It would now be interesting to explore 
such an adaptive response experimentally. For example, such an oscillatory adaptation was 



recently observed in cAMP concentration in Dictyostelium cells [21 1. 

Generally, cells undergo autocatalytic reactions to replicate themselves, where complex 
autocatalytic reactions are often slower than simple catalytic ones. The results of this study 
suggest that adaptation obeying Weber's law is generally observed in such systems with 
slow autocatalytic reactions. For example, consider a sequential reaction process for the 
synthesis of a biopolymer with catalytic activity. With an increase in the sequence length 
of the polymer, more time would generally be required for its synthesis and degradation; 
whereas, polymers with a longer sequence have the ability to catalyze these chain reactions. 
This process is in agreement with the sequential autocatalytic model studied here. 

The present model does not show FCD in the strong sense [12J]. However, the long-term 
relaxation process depends only on the fold change and not on the absolute value of the 
external parameter S, under the adiabatic condition. 
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Adaptation or habituation refers to a general property of a biological system to exhibit 
homeostasis, whereas Weber's law is applicable to a system exhibiting a sensory response 
over a wide range of environmental conditions. The variation in an external parameter S 
will often increase with an increase in its value. Then, if the response were just proportional 
to the difference in the external parameter before and after the input, it would be too 
sensitive to large S and would not be able to generate a reliable response under external 
noise. Therefore, a response to fold change in accordance with Weber's law is applicable to 
sensory systems in general. 

In the present study, Weber's law was usually obeyed in an autocatalytic reaction system 
with timescale differences, without any special design or tuning of the parameters. Because 
of its simplicity, the present mechanism is expected to have a wide range of applications. 
It may also offer fresh perspectives on Weber's law in general, including psychological and 
neural perspectives. Indeed, our proposed model can be regarded simply as a system with 
self-positive feedback and is not necessarily restricted to chemical reactions. 

The authors would like to thank K. Kamino, S. Sawai, and K. Fujimoto for useful dis- 
cussions. M.I. was partially supported by the JSPS. 
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